

IONOSPHERIC RESEARCH FROM SPACE VEHICLES 


*•3 '212 80 


R. E. BOURDEAU 
Space Sciences Division 

NASA-Goddard Space Flight Center, Greenbelt, Maryland, U.S.A. 


(Received January 12, 1963) 

1. Introduction 

Previous to the advent of artificial satellites, our knowledge of the ionosphere was 
limited principally to the interpretation of data from the classical ground-based iono- 
sonde and a few rocket measurements of ionizing radiation, electron density and ion 
composition, all applicable to altitudes below the F2 peak. The early rocket results of 
Seddon and Jackson, which have been summarized by Ratcuffe and Weekes 
(1960), show that the lower ionosphere is characterized by an electron density which 
increases monotonically with altitude and thus that the term “layer” is an incorrect 
nomenclature for the D, E, and F regions. In more recent years, a wealth of data has 
been gathered on the temporal, spatial and energy distribution of the charged particles 
which for the ionosphere. This has been accomplished at previously relatively un- 
explored altitudes well above the F2 peak. Concurrently but unfortunately not 
simultaneously, measurements have been made of the time dependence of corpuscular, 
ultraviolet and X-radiation. The most serious observational gap is in the structure of 
the neutral atmosphere where the most significant contribution has been the computa- 
tion of the neutral density parameter from satellite drag observations. 

This paper compares the available interdisciplinary experimental data with theo- 
retical models of the D, E, and F regions as well as the upper ionosphere, eventually 
selecting those models which best fit the spaceflight observations. Because a significant 
amount of low-energy charged particle data were obtained by direct sampling 
techniques and because this methodology is relatively new to ionospheric research, the 
paper contains brief discussions of the validity of some of these results as they are 
presented. Although the theoretical models of all ionospheric regions have been 
enhanced significantly by increased spaceflights, refinement and changes in these 
models await the next two major steps - the launching of satellites and rocket probes 
which are truly geophysical in nature and the correlation of such interdisciplinary 
measurements with data resulting from recent breakthroughs in ground-based 
observational methods. 

2. The D Region Under Quiet Solar Conditions 

The D region, situated between approximately 50 and 85 kilometers, is the lowest 
! ionospheric subdivision where a significant number of free electrons are found. This 
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altitude interval is difficult to treat theoretically and experimentally because of the low 
charged particle concentrations, because the relatively high gas densities result in high 
electron 'collision frequencies, and because of the high probability of negative ion 
formation. 

Spaceflight radiation observations made during the absence of solar flares show 
that the three most probable ionizing agents are cosmic, Lyman a (1215.7 A) and 
X-radiations (2-8 A). The most detailed theoretical analysis of D region formation 
has been performed by Nicolet and Aikin (1960) using an expansion of the following 
equation of ionization: 

q (oc d + boti) n+n e , (1) 

where q is the equivalent electron production rate, oca is the loss coefficient for re- 
combination of positive molecular ions with electrons, a* the loss coefficient for re- 
combination of positive with negative ions, b is the ratio of negative ions to electrons, 
and n e and n+ are the electron and positive ion concentrations, respectively. The cosmic 
and X-radiations act on molecular nitrogen and oxygen, the major neutral consti- 
tuents, whereas Lyman a radiation acts on nitric oxide, a trace constituent. In their 
expansion of eq. (1), Nicolet and Alkin produced a summation of the separate effects 
of the three ionizing sources and concluded that: 

(a) cosmic radiation is the most important ionizing agent at altitudes up to 70 km ; 

(b) in the 70-85 km region, assuming an N0+ recombination coefficient of 
3 X 10 9 cm 3 sec -1 , a nitric oxide abundance of only 10 10 of the total neutral 
concentration is required to make X-rays (2-8 A) unimportant to the formation of the 
D region under quiet solar conditions ; 

(c) negative ions are unimportant above 70 km ; 

(d) the ionization at altitudes between 85 and 100 km (base of the E region) is the 
result of X-rays (30-100 A) and ultraviolet radiation (Lyman /? and the Lyman 
continuum). 

It is known principally from measurements made on the Greb satellite (Kreplin, 
Chubb and Friedman, 1962) and on the Orbiting Solar Observatory (Lindsay, 1962) 
that the Lyman a energy flux (3-6 erg cm -2 sec x ) is relatively constant with solar 
condition. The X-ray energy flux, however, is extremely variable. Consequently, it is 
possible that Lyman a and X-radiations alternate as the predominant ionizing agents 
of the normal upper D region. Poppoff and Whitten (1962) have challenged the 
Lyman a hypothesis in attempting to account for the 70-85 km ionization entirely by 
X-radiation. However, they used energy flux values, which according to a rcente 
review (Friedman, 1962), exceeds energies representative of a quiet sun even at the 
maximum of the solar cycle. 

In this section it will be shown that to achieve agreement with spaceflight n+ ob- 
servations taken during quiet midsolar cycle conditions, it is necessary to invoke the 
Lyman a hypothesis. Additionally it will be shown that to obtain consistency with 
ground-based measurements of electron density for similar conditions, existing D 
region formation theories must be considered incomplete. Specifically, itappears that 
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there is a tendency in most of the current models to underestimate the number of 
electrons lost in the various processes of negative ion formation. * 

It has been inferred in the previous discussion that care must be taken in comparing 
experimental results with theoretical models for the same solar condition. The extreme 
sensitivity of the X-ray energy flux to solar activity is demonstrated in Figure 1. It is 



Fig. 1. The variability of solar X-ray emissions. Curves A, C and E are from Friedman (1962) 
Curves B and D are from Pounds and Willmore (1962). 

seen by a comparison of curves A and C (Friedman, 1962) that even in the absence of 
flares the intensity in the important wavelength region (2-8 A) varies by two to three 
orders of magnitude during a complete solar cycle. The recent data (curve B) from 
the Ariel satellite (Pounds and Willmore, 1962) correspond temporally to values 
used by Nicolet and Aikin (1960) in computing their “quiet sun” model. The Ariel 
data also correspond most closely to the epoch of the solar cycle under which the 
few available D region charged particle density profiles have been obtained. 

In Figure 2 are presented two theoretical electron density profiles (Aikin, 1962a) 
which are intended by their contrast to illustrate the enhancement of n e that should 
occur in the 70-85 km region if the NO concentration is 10~ 10 of the total concentration 
and if the negative ion abundance is that computed with the assumption that O 2 is the 
principal negative ion. The alternative theoretical profiles were computed by using 
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X-ray energy fluxes approximately representative of curve B in Figure 1. The principal 
change from the original Nicolet-Aikin model is that the values of oca for Oj andNO + 
were increased by a factor of ten. The Oj coefficient is increased in the light of more 
recent laboratory investigations (Kasner, Rogers and Biondi, 1961) whereas in a 
succeeding section it is shown that an increase of ad for NO + from 3 x 10 -9 to 



Fig. 2. Electron density profile of the quiet daytime D region at mid-latitudes for the period 
between solar maximum and solar minimum. 


3 x 10~ 8 cm 3 sec -1 is more consistent with an effective recombination coefficient 
determined for the E region from diurnal variation observations. 

To date, a reliable spaceflight method for measuring the low electron densities 
characteristic of the normal D region has not yet been developed. USSR data from a 
rocket-borne impedance probe for altitudes as low as 75 km have been published 
(Krasnuskin and Kolesnikov, 1962) but too recently for critical review. The lag in 
the development of rocket-borne electron density experiments is due primarily to the 
extremely low electron to neutral gas density ratio and the associated high collision 
frequencies. Because of this lag, it is necessary to resort primarily to electron density 
profiles obtained by ground-based methods. 

In a recent review, Ratcliffe and Weekes (1960) conclude that D region electron 
density profiles obtained by ground-based methods should be treated with caution. 
According to this review, the most reliable of these profiles was obtained by Gardner 
and Pawsey (1953) using observations of weak backscatter echoes. The experimental 
data were obtained at a time when the values of the X-ray energy flux most likely were 
intermediate between those represented by curves A and B in Figure 1. The experi- 
mental data are plotted in Figure 2 for comparison with the alternative theoretical 
models. The valley of ionization which shows a minimum at about 78 km is typical of 
most of the few electron density profiles reported by other observers using ground- 
based techniques. Belrose and Burke (1961), for example, report a valley with the 
minimum located at 70 km, also using a backscatter method. There is a large difference 
between both theoretical ionization estimates and the experimental data at altitudes 
below 70 km. It has been concluded by others (Belrose and Burke, 1961) that 
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observational data by the backscatter method must be treated with caution in the 
lower D region because of simplifying assumptions in the Afpleton-Hartree 
formula used to derive the experimental result. If the experimental data below 70 km 
are assumed correct, the low electron densities could be explained either by a theo- 
retical overestimate of the effective cosmic ray electron production rate or an under- 
estimate of the processes leading to the formation of negative ions or a combination 
of the two. 

It is seen by a comparison of the experimental data with the two alternative 
theoretical models in Figure 2 that all of the ionization in the 70-85 km region could 
be explained from the X-ray effect alone as proposed by Poppoff and Whitten with- 
out resort to the Lyman a hypothesis of Nicolet and Aikin. However, a conclusion 
based on these few electron density data that the Lyman a radiation has a negligible 
effect on the formation of the D region would not be consistent with spaceflight 
observations of the positive ion characteristics, a discussion of which follows. 

The first spaceflight measurements of D region ion parameters (Bourdeau, 
Whipple and Clark, 1959) were obtained by use of a rocket version of a Gerdien 
condenser similar to that flown extensively on aircraft and balloons by workers in 
atmospheric electricity. This experiment measures ion conductivity (A+). In the al- 
titude region where the number of free electrons is negligibly small, A+ can be estimated 
theoretically from 

l+ = ek+(q+/* T )W f (2) 

where e is the elementary charge, k+ is the ionic mobility, q~ is the rate of ion-pair 
production and ay is Thomson’s ion-ion volume recombination coefficient. The 
altitude dependence of k+ and olt are estimated from laboratory experiments. In the 
reported results, the altitude at which <x T can be extrapolated with confidence was 
overestimated. However, at altitudes up to 50 km where the use of ay is appropriate, 
the experimental values of /+ are in good agreement with theoretical estimates based 
on Equation (2) in which q+ is estimated from spaceflight observations of cosmic ray 
intensities. The experimental rocket values of 2+ also are consistent with those obtained 
on balloons by other investigators. The importance of these agreements is that it lends 
some confidence to the theoretical conclusion (Bourdeau et al ., 1959; Ichimiya, 
Tokayama and Aono, 1959) that possible errors in the rocket experimental data 
associated with shock waves (thermal ionization, adiabatic compression etc.) are 
small, even at the relatively high D region gas pressures. 

Positive ion densities now have been reported by Smith (1961a) and Aono, Hirao 
and Miyazaki (1961), each using a different method of computing n+ from the ob- 
served current to an exposed electrode insulated from the rocket body. Their published 
values are compared in Figure 3 with a theoretical n+ model (Aikin, 1962a). Estimates 
of negative ion densities based principally on attachment and photodetachment 
processes are represented by the difference in the theoretical n + and n e curves. Both 
sets of experimental data were obtained at a time when the values of the X-ray energy 
flux were most likely slightly larger than those represented by curve B in Figure 1. 
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It is seen from Figure 3 that the shape of the ion density profile experimentally 
obtained by Smith is consistent with the shape of the theoretical n+ curve and thus 
consistent with a hypothesis of three different ionization sources between 50 and 
100 km, each predominating in a separate altitude region. However, it is seen also that 
the experimental values of n+ are generally much higher than the theoretical values. 
The general excess of the experimental positive ion densities over the theoretical 
estimates may be interpreted as an altitude invariant error either in the experimental 



Fig. 3. Ion density profile of the quiet daytime D region at mid-latitudes for the period between 

solar maximum and solar minimum. 


result or in one of the assumed theoretical parameters. It now will be shown that 
experimental uncertainties could account for the difference between the theoretical 
and experimental n+ curves in Figure 3. 

The results reported by Smith (1961a) were obtained by the use of a standard 
asymmetric Langmuir probe in which the active electrode is the nose tip of the rocket. 
In the D region this Langmuir probe functions as a conductivity meter, the principle of 
operation being analogous to that of the exploratory Gerdien condenser experiment 
(Bourdeau et al 1959). Positive ion current flowing from the ionosphere to the nose 
tip is measured as a function of a variable voltage applied between this electrode and 
the main rocket body. Values of A+ are computed by electrostatic theory from the slope 
of the volt-ampere curve when the active electrode is negative with respect to the 
ambient plasma. Smith computed positive ion densities from the observed A+ accord- 
ing to the following simplified form of the pertinent electrostatic equation: 


n+ = (A+m+v+)/e 2 , (3) 

where m+ is the mean ionic mass and v+ is the ion collision frequency. Laboratory 
investigations (Loeb, 1955) indicate that the simplifying assumptions leading to Equa- 
tion (3) could result in uncertainties of a factor of two in the n+ determination. Smith 
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used values for v+ contained in the 1959 ARDC model atmosphere which are probably 
uncertain by another factor of two (Aikin, private communication). Consequently, 
the values of n+ computed by Smith could be uncertain by a factor of four just from 
the geophysical interpretation of the measured conductivity alone. This uncertainty is 
greater at altitudes above 80 km because the methods used to obtain A+ from the 
observed positive ion current are no longer valid when the ionic mean free path 
becomes large compared to the dimensions of the active electrode (Bourdeau et al , 
1959). 

Other possible errors in the measurement of D region ion characteristics are 
associated with aerodynamic effects (discussed previously) and the photocurrent 
effect. The latter effect represents the predominant criticism of ionospheric direct 
measurements techniques. For this reason, it is necessary to present the following 
arguments for a neglible photocurrent effect on the experimental n+ results shown in 
Figure 3. 

Because solar radiation can induce electron emission from rocket surfaces, an 
artificial electron current will flow from this electrode into the surrounding medium. 
This unwanted photocurrent is of the same polarity as the desired positive ion current. 
Values for the photocurrent density ( j p ) of 2.3 x KF 9 ampcm -2 at an altitude of 
200 km have been measured on a rocket (Hinteregger, Damon and Hall, 1959) and 
of 5 X 10 -9 ampcm -2 in the upper ionosphere on the Explorer VIII satellite 
(Bourdeau, Donley, Serbu and Whipple, 1961). The approximate agreement of the 
j p values measured at such radically different altitudes should not be surprising 
because the most responsible solar radiation, occurring at wavelengths between 2000 
and 3000 A (Hinteregger et al , 1959; Ichimiya et al , 1959), is not significantly 
absorbed even in the D region. 

The values of the perturbing photocurrent density listed above can be shown to be 
some two or three orders of magnitude larger than the ambient random ion current 
density flowing in a D region undisturbed by the presence of a rocket. However, the 
actual effect is dependent on the ratio a+ja p where a+ is the effective ion collection 
area and a p is that part of the active electrode’s surface area which is exposed to the 
sun. This ratio is extremely difficult to obtain theoretically for the D region although 
orders of magnitude at large negative electrode potentials are indicated (Ichimiya 
et al , 1959). It is possible by this mechanism then, to overcompensate in magnitude 
for the photocurrent effect. Confidence in this conclusion can be acquired by examina- 
tion of observed D region volt-ampere curves. At negative electrode potentials, the 
current measured at the electrode is observed to be an order of magnitude larger than what 
would be expected from the product of j p and a p . Additionally even if the amplitude 
of the photocurrent was large, /+ is obtained from the slope of the volt-ampere curve 
at negative electrode potentials. This is a regime where the photocurrent is invariant 
with electrode potential. Hence it should not affect the measurement of the slope and 
thus the n+ determination. 

The n+ data by Aono et al (1 961 ) shown in Figure 3 were obtained by measurement 
of the positive ion current to an effectively spherical electrode which was extended 
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along the major rocket axis. The term “effectively spherical” is used because although 
the electrode actually consists of two crossed rings, they are so arranged that the 
analysis for ions is as though the electrode were a solid sphere. The virtue of the 
crossed rings is that they present a much smaller target to solar radiation and con- 
sequently the experiment is even less sensitive to the photocurrent effect than the 
device used by Smith. 

A form of the conventional Langmuir probe theory applied to ions is used by the 
Japanese in computing n+ from the experimental volt-ampere curves. The general 
agreement of their published n+ values with an independent simultaneous measurement 
of n e in the E region speaks well for the validity of their application of this theory at 
least for the simple case of no collisions in the plasma sheath surrounding the probe. 
However, their published data show slightly larger electron than ion concentrations 
in some portions of the E region which are in turn indicative of second-order errors 
even for the simple case. The validity of Langmuir probe theory becomes more un- 
certain at D region altitudes where multiple collisions within the sheath complicate 
the application of the conventional Langmuir probe theory. Thus, for different reasons 
associated with the interpretation of the observed volt-ampere curves the largest un- 
certainty in both sets of experimental data shown in Figure 3 can be expected between 
80 and 90 km. It would not be surprising, therefore, if the same factor of 4 in the 
difference between the experimental data in the 80-85 km region were obtained even 
under identical solar conditions. 

It generally is concluded by this author that the uncertainties in the D region 
experimental n+ values could be large enough to adjust to the theoretical n+ estimates 
of Aikin (1962a). On the other hand, it is unlikely that the experimental inadequacies 
are sufficient to adjust rocket observations of n+ to the experimental electron density 
values shown in Figure 2, thus suggesting a much higher negative ion abundance than 
estimated theoretically by Nicolet and Aikin (1960) and by Whitten and Poppoff 
(1961). The high negative ion abundance is more in accordance with the theoretical 
model of Moler (1960). The suggestion of high negative ion densities also is consistent 
with values computed from the original experimental conductivity data by Whipple 
(1960) who suggests that dust particles believed responsible for noctilucent cloud 
formation at high latitudes could provide an important recombination surface for 
charged particles in the D region. The simultaneous measurements of electron and 
positive ion densities made by Aono et al. (1961) show somewhat larger n+ values in 
the 95 km region again suggesting a large negative ion abundance at the lower D 
region altitudes where their formation is even more likely. 

The few data that are available support a tentative conclusion of ionization of the 
lower and upper D regions by cosmic and Lyman a radiation, respectively, during 
conditions of a quiet sun at the middle of the solar cycle. A firm conclusion awaits 
more accurate measurements of electron and ion densities obtained simultaneously 
with observations of the pertinent energy fluxes of the ionizing sources. Just as 
important is the need for improved knowledge of the pertinent reaction rates and of 
the number density of the trace neutral constituents. 
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3. The D Region During Disturbed Solar Conditions 

During periods of high solar activity, the D region is characterized by enhanced 
ionization with associated electromagnetic wave attenuation strong enough to produce 
radio blackouts. There are many phenomena associated with solar flares which in- 
crease the normal D region electron densities by up to two orders of magnitude. A 
comparison of the amount of enhanced ionization for different types of solar flare 
events is made with the Nicolet-Aikin quiet sun model (curve A) in Figure 4. 



Simultaneously with the appearance of the flare, radio absorption is observed in 
the D region on the sunlit side of the earth for periods lasting up to approximately one 
hour. This particular phenomenon is called a Sudden Ionospheric Disturbance. 
Chubb, Friedman and Kreplin (1960) made rocket flights into such events and 
observed abnormally high X-ray fluxes at extremely low D region altitudes. These 
observations have been used in Equation (1) by Nicolet and Aikin (1960) to calculate 
an S.I.D. electron density profile (curve B of Figure 4). Here, the influences of cosmic 
and Lyman a radiation become minor, the secondary layer in the 70-85 km region dis- 
appears, and the profile is characterized by overall enhancement and a monotonically 
increasing electron density. The S.I.D. profile does not include the contribution of 
radiation responsible for the formation of the E region. Observations from the Greb 
satellite show that whenever the X-ray flux at wavelengths shorter than 8A exceeds a 
critical value of 2 x 10~ 3 erg cm -2 sec -1 , radio fadeouts at short wavelengths and 
other sudden ionospheric disturbances occur (Friedman, 1962). It is seen from curves 
D and E in Figure 1 that this high flux value generally would be characteristic of flares 
stronger than importance 2. 

A second type of ionospheric storm predominates at night and is associated with 
active aurorae and magnetic disturbances. It generally is accepted that the ionizing 
agents are energetic particles comprising corpuscular emission from the sun, a belief 
founded on the observation that D region absorption occurs some 21 hours after the 
appearance of a flare, an interval corresponding to the sun-earth transit time for these 
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particles. During the storm, D region electron densities increase to values high enough 
that echoes are observed on ionosondes. ‘‘Layers” appear at 90 km during weak geo- 
magnetic activity and as low as 70 km for the more active events. 

A third type of disturbance occurs only at auroral latitudes. These phenonema, 
called Polar Cap Absorption (PCA) events, are produced by energetic protons 
emitted from the sun during certain flares. Here echoes are observed from ionosondes 
at altitudes as low as 60 km. The phenomenon has been the subject of considerable 
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Fig. 5. Rocket measurement (Kane, 1 960) of the altitude dependence of electron collision frequency. 


study during the past few years. Most recently, Maehlum and O’Brien (1962) have 
proposed a semi-quantitative time history of the D region electron density profile 
during a PCA by inserting proton fluxes observed on the Injun satellite as a source 
function in the conventional equation of ionization. One of these profiles (curve D of 
Figure 4) coincides with the peak of 27.6 Me cosmic radio noise absorption (17 db) 
observed at College, Alaska. Like the Nicolet-Aikin model for an S.I.D., curve D 
does not contain the contribution due to solar radiation. Inclusion of solar radiation 
should modify the result in the base of the E region (above 80 km). Jackson and 
Kane (1959) used a rocket-borne radio propagation experiment to measure a D region 
electron density profile (curve C of Figure 4) during a less active phase of a PCA event 
(3 db absorption at 30 Me). This profile compares favorably in shape but is generally 
less enhanced than that proposed by Maehlum and O’Brien for an interval when the 
observed radio absorption also was 3 db. If, as the evidence presented in the previous 
section indicates, there is a higher D region negative ion abundance than is generally 




IONOSPHERIC RESEARCH FROM SPACE VEHICLES 


693 


accepted, then the three dashed curves in Figure 4 could represent overestimates of the 
electron density. 

The strong absorption of radio waves in the D region is made possible by the high 
electron collision frequencies characteristic of this altitude interval. Kane, (1959, 
1960) has derived an electron collision frequency model from measurements of the 
difference in the absorption of ordinary and extraordinary radio propagation modes. 
His data which represent the best available rocket measurements of the altitude 
dependence of this important D region parameter are presented in Figure 5. They 
were obtained on two separate rocket flights flown during auroral absorption events. 

4. The Ion Content of the Lower Ionosphere 

To a high degree of probability, the predominant ionizing sources of the E and F 
regions for quiet solar conditions at mid-latitudes are solar ultraviolet and X-ray 
radiations. It is possible to compute the rate at which different ion species are formed 
from a knowledge of the altitude dependence of the number of photons incident at 
each wavelength, of the densities of the individual neutral constituents and of the 
absorption and ionization cross sections of these constituents. There are available 
some altitude profiles of the photon fluxes as a function of wavelength but only for 
limited conditions. Early results from rocket-borne neutral gas spectrometers re- 
ported to date are extremely controversial. The validity of these results is questioned 
justifiably because recombination effects within the instruments distort the gas under 
study from its ambient condition. Thus most model atmospheres do not depend on 
experimental observations but rather are derived from the hydrostatic laws using an 
assumption for the critical altitude of diffusive separation. Our knowledge of the 
important cross sections also must be considered incomplete. 

The most recent estimate of the altitude dependence of the rate at which various 
ions are produced was made by Watanabe and Hinteregger (1962) using solar 
radiation data obtained by Hinteregger under quiet solar conditions at mid- 
latitudes. The rates for the ions formed in greatest abundance (Ng, 0 % 0 + ) are illus- 
trated in the left-hand side of Figure 6. Watanabe and Hinteregger hasten to point 
out that these curves must be considered suggestive rather than quantitative princi- 
pally because of the uncertainty in our knowledge of the neutral atmospheric composi- 
tion. Some adjustment of these rates already are indicated from very recent neutral gas 
spectrometer results (Schaeffer, 1963). These results were obtained by use of a spec- 
trometer which overcomes the inadequacies of the early devices by ionizing the gas 
under study before its entry into the instrument’s analyzer section. The results show a 
larger ratio of atomic to molecular oxygen than that used to obtain the ion production 
rates shown in Figure 6. 

Some indication of the solar radiation responsible for the production of various 
ionospheric regions can be made by comparing the inferred production rates of the 
individual ionic species with rocket-borne ion spectrometer results. Early rocket 
flights (Johnson, Meadows and Holmes, 1958) of a radio-frequency mass spectro- 
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meter made in the auroral zone showed that the diatomic ions, Oj and NO + , pre- 
dominate below 200 km and that 0 + is the principal ion found above this altitude. 
These results have been confirmed by Istomin (1960) and by Taylor and Rrinton 
(1961). The rocket data of Taylor and Brinton shown in the right hand side of 
Figure 6 are selected for comparison with the production rates because they were 
obtained for latitude and temporal conditions most closely representing those under 
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Fig. 6. Comparison of the production rate and the actual abundance of ions in the lower ionosphere. 


which the photon fluxes used in estimating the photoionization rates were obtained. 
It is seen in the comparison that although the diatomic nitrogen ions are expected to be 
produced in great quantities, the spectrometer results show that they exist only as 
extremely minor charged constituents. The most likely reason is that Nj ions dis- 
sociatively recombine (Nj + e -> N' + N") at a very rapid rate as demonstrated by 
Faire and Champion (1959). 

It was inferred in the previous section that there is general acceptance that X- 
radiation at wavelengths longer than 10 A and ultraviolet radiation (Lyman ft and the 
Lyman continuum) are responsible for the base of the E region (85-100 km). There 
remain two schools of thought for the formation of the remainder of the E region: 
(a) general ionization of air by soft X-rays (Vegard, 1958; Hulburt, 1958; Bates, 
1956; Friedman, 1959) and (b) photoionization of O 2 by solar radiation in the wave- 
length region between 800 and 1026 A (Wulf and Deming, 1938; Nicolet, 1945; 
Watanabe, Marmo and Pressman, 1945). On the basis of more accurate solar radia- 
tion data, Watanabe and Hinteregger (1962) strongly support ionization of the 
100-120 km region by direct production of Oj from ultraviolet radiation at wave- 
lengths between 911-1027 A. Part of this conclusion is based on the agreement 
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between experimentally observed E region electron densities and those computed 
from the conventional ionization equation under quasi-equilibrium conditions: 

Tie = ( q/oid ) 1/2 , (4) 

where q is the photoionization rate for Oj computed from the observed ultraviolet 
photon fluxes. 

However, there are two sets of recent experimental evidence which indicate that 
Watanabe and Hinteregger overestimated the electron density produced in the 
100-120 km region by ultraviolet radiation. Firstly, the recent rocket neutral gas 
spectrometer results indicate that the O 2 concentration and thus the Oj production 
rate have been overestimated. Secondly, Watanabe and Hinteregger used a value 
for the Oj recombinance coefficient taken from Nicolet and Aikin (1960) of 3 x 10 -8 
cm 3 sec -1 which as discussed in the previous section appears from more recent 
laboratory work to be an order of magnitude too low. 

A third set of experimental data which also tends to favor X-radiation as the more 
important agent for E region formation involves the rocket ion spectrometer data 
shown in Figure 6. The diatomic ion NO+ is observed to be predominant at E region 
altitudes. Its existence can best be explained (Aikin, 1962b) by ion-atom interchange 
involving the 0 + ion (O 4 + N 2 NO + + N). The ultraviolet radiation (91 1-1027 A) 
observed to penetrate below 120 km (Hinteregger and Watanabe, 1962) lies in a 
range above the ionization threshold for 0 + . Consequently to be consistent with the 
rocket ion spectrometer observation that NO + is the predominant E region ion, the 
ultraviolet hypothesis requires a reaction leading to NO + formation involving Og, a 
reaction not considered by Hinteregger and Watanabe. Parenthetically, it is difficult 
to dispute the correctness of the observed ratio NO + /C >2 from an experimental point 
of view since the principal corrections which must be made to the measured ion 
spectrometer target currents involve the effects of vehicle motion and of electric fields 
due to the attracting vehicle potential. Both of these effects are expected to be nearly 
the same for ions of approximately the same mass. Hence, the measurement of their 
relative abundance should be reasonably accurate. As Hinteregger and Watanabe 
(1962) demonstrate, the ion spectrometer measurement of the increasing importance 
of 0 + with altitude leading to its predominance above 200 km is explained adequately 
by the simultaneously increasing role of extreme ultraviolet radiation in the wave- 
length region between 280 and 911 A. 

5. The Diurnal and Anomalous Behavior of the E Region 

The altitudes generally assigned to the E region lie between approximately 85 and 
140 km. The most accurate rocket method of daytime electron density determination 
for this altitude region remains the Seddon two-frequency radio-propagation experi- 
ment. Results from this experiment early demonstrated that the lower ionosphere is 
not characterized by discrete layers but rather that the electron density increases mono- 
tonically with altitude and thus that the F region is a continuation of the upper E 
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region. Until recently, there have been no detailed observations of nighttime E region 
ionization because neither this rocket-borne experiment nor ground-based ionosondes 
are sensitive enough. Now the introduction of low energy plasma probes into space- 
flight studies of the ionosphere permits our first insight into the diurnal characteristics 
of the E region. 



Fig. 7. Rocket measurements of the diurnal and anomalous behavior of the E region. 


Electron density values reported from one daytime and two nighttime rocket 
launchings of an asymmetric Langmuir probe (Smith, 1961a; Smith, 1961b) are 
presented in Figure 7. It is seen that the electron density in the lower portion of the 
nighttime E region is found to be 3 X 10 3 cm -3 at approximately 22h LMT and 
1 x 10 3 cm 3 near sunrise, both about two orders of magnitude less than typical 
daytime values. Smith (1961b) has used this observed diurnal variation to compute an 
effective E region recombinance coefficient of 2 x 10 -8 cm 3 sec -1 , which is perhaps 
mostly representative of the observed predominant ion, NO+. In the view of this 
author, this value should be considered approximate, firstly because of possible 
second-order errors in computing absolute values of n e from the observed volt- 
ampere curves and secondly because it is assumed that there is no nighttime ionization 
source. It is tempting from the observation of deep ionization valleys characteristic of 
the upper nighttime E region profiles shown in Figure 7 to suggest that the ledges in 
the lower E region result from observed scattering (Donahue, 1962) of ultraviolet 
radiation in the night sky. 

The other significant data contained in Figure 7 are associated with the measure- 
ment of the ionization characteristics of the sporadic-E layer (£ s ), which on an iono- 
sonde is characterized by echo reflection from a constant virtual height. The statistical 
occurrence of the E s phenomenon from ionosonde data has been the subject of 
considerable study, yet the causes of sporadic-E layers remain topics for speculation. 
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Early rocket observations of daytime E s ionization were obtained by Seddon and 
Jackson (1958). 

Two types of E s ionization, differing in horizontal dimensions are contained in 
Figure 7. It is seen that the electron density observed in the sporadic-E layer (102 km) 
on the 17 August flight represents an enhancement of about a factor of four above the 
average lower E region value. The thickness of the layer measured at half the peak 
electron density was 450 meters. The layer was detected at precisely the same altitude 
on both the ascent and descent portions of the rocket trajectory, indicating a horizon- 
tal dimension greater than 72 km. Similar characteristics were observed at an altitude 
of 112 km on the 27 October flight. It has been suggested (Whitehead, 1960) that 
E s ionization with these characteristics results from re-distribution of electrons by 
wind shear in the presence of a geomagnetic field, rather than local changes in electron 
production or loss. 

On the 27 October flight, another type of E s ionization was observed. One such 
electron density peak is shown at an altitude of approximately 108 km in Figure 7. A 
larger number of such peaks than indicated in Figure 7 were observed with no correla- 
tion between the ascent and descent portions of the trajectory. The small horizontal 
dimensions lead Smith (1961b) to support current theories that some types of E s 
ionization result from meteor trails. 

It can be said that there is insufficient evidence to support conclusively any current 
theories to account for E s ionization. Such evidence will come from planned correlative 
experiments performed simultaneously with charged particle density measurements. 
The rocket data are, however, serving to define the dimensions and ionization char- 
acteristics of sporadic-E layers. One of the most detailed altitude definition of the E s 
ionization reported to this author was obtained by Boyd (private communication) 
from a daytime ion density profile on December 6, 1961. His observation of a depth of 
just under one kilometer and an ionization enhancement of about a factor of 3 is 
consistent with the results shown in Figure 7. 

In addition to the measurement of electron density the asymmetric Langmuir 
probe is designed to provide data on the electron temperature, T e . For the same rocket 
flight from which the daytime electron density profile illustrated in Figure 7 was ob- 
tained, Smith (1961a) reported electron temperature values of 1400° K with little 
altitude variation between 100 km and the apogee altitude of 155 km. The disagreement 
of these T e values with current hypotheses for this altitude region and with the results 
of other rocket experiments made under quiet daytime conditions justifiably has 
raised some questions as to the validity of the electron density results. There are two 
possible explanations for the differences. The first, geophysical in nature, is that this 
particular rocket flight took place within twenty-four hours of a geomagnetic 
disturbance. Current theoretical concepts of expected differences between the electron 
and neutral gas temperatures do not include the extreme variability of X-ray fluxes 
which perhaps could account for the high electron temperature. On the other hand it is 
just as likely, from this author’s point of view, that the high electron temperatures 
result from simplifying assumptions used in deriving T e from the observed volt- 
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ampere curves. Specifically, n e and T e both are derived from the volt-ampere curve 
after unwanted photocurrent and positive ion currents are subtracted out. T e is 
computed from the slope of the corrected curve at negative collector potentials and is 
sensitive to the accuracy of the correction. The electron density, on the other hand, is 
computed from the current observed at a distinct discontinuity in the curve where the 
collector is at plasma potential. At this point the electron current is more than an 
order of magnitude higher than the unwanted current and hence the derivation of n e is 
insensitive to the assumed corrections. Thus the errors in deriving n e are small enough 
to have little effect on the important conclusions from Figure 7 relative to (a) the 
magnitude of the diurnal variation of n e in the lower E region, ( b ) the observation of a 
deep ionization valley in the upper E region at night, and (c) the major characteristics 
of sporadic-E ionization. 


6. The Formation of the F2 Peak 

The altitudes generally assigned to the F region lie between 140 km and the altitude at 
which 0 + ceases to be the predominant ion. This region contains the altitude (F2 peak) 
of the maximum electron density found in the ionosphere. It was not until very recent 
years that significant experimental data on the characteristics of the upper F region 
were obtained. Local charged particle density measurements now have been made on 
satellites using ion traps (Krassovsky, 1959; Bourdeau, 1961) and rf impedance 
probes (Bourdeau and Bauer, 1962; Ulwick and Pfister, 1962; Sayers, Rothwell 
and Wager, 1962). 

More meaningful to the physics of the F region have been experimental measure- 
ments of the charged particle density distribution with altitude where latitudinal and 
temporal variations can be neglected. Such profiles have been obtained from ground- 
based experiments (Van Zandt and Bowles, 1960), rocket-borne radio propagation 
experiments (Berning, 1960; Jackson and Bauer, 1961 ; Knecht and Russell, 1962) 
and ion traps (Hanson and McKibbin, 1961 ; Hanson, 1962a). These investigators all 
have concluded from their observations of an electron density distribution with a 
practically constant logarithmic slope taken over a few hundred kilometers that the 
upper ionosphere can be represented as an isothermal medium in diffusive equilibrium. 

The most favored theory of the formation of the F2 peak at mid-latitudes under 
quiet solar conditions can be illustrated by comparing these charged particle density 
profiles with the altitude dependence of electron production rate inferred by Watanabe 
and Hinteregger (1962). Such a comparison is made in Figure 8, using the electron 
density profile obtained by Jackson and Bauer (1961) with a rocket-borne cw radio 
propagation experiment. It has been recognized that, as illustrated, the altitude of the 
F2 peak lies considerably above that of the maximum electron production rate. This is 
best explained from the following form of the continuity equation corresponding to 
quasi-equilibrium : 


d 

q = pn e + (n e W d ) , 
dr 


( 5 ) 
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where fi represents an attachment - like loss coefficient, z is the altitude and W d serves 
to define changes in electron density by vertical motion or diffusion. In the region up 
to the F2 peak, ftn e predominates over the diffusion term and because it decreases more 
rapidly with altitude than the electron production rate, q , n e increases with altitude up 
to its maximum value. Near the F2 peak, the diffusion term becomes larger than j3n e 




Fig. 8. Comparison of the altitude dependence of electron production rate and electron density, 
illustrating the formation of the F2 peak. 

and the electron density begins to decrease. The diffusion mechanism is caused by 
gravitational forces acting upon the ions which, by coulomb attraction, cause the 
electrons to diffuse downward. At altitudes well above the F2 peak where q and fin e in 
Equation (5) can be considered negligible, the electron density distribution is given by 
the diffusive equilibrium equation: 

d(ln n e ) = m+g 

dz ~ kiTe+ Ti)' U 

where g is the acceleration of gravity, k is the Boltzmann constant, and T e and 7* are 
the electron and ion temperatures. This most favored model of the formation of the 
F2 peak perhaps requires some modification near the geomagnetic equator where 
because of the existence of a horizontal magnetic field, diffusion cannot readily 
produce vertical movements. 

Many empirical attempts have been made to solve the continuity equation, 
Equation (5), by making assumptions about the nature of the vertical movement, the 
production term, q , and the loss coefficient, /?, and comparing the calculated results 
with experimental ones. Principally because of the numerous variables involved, these 
models vary considerably - from attempts to explain the F2 peak without diffusion 
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(Sagalyn and Smiddy, 1963) to agreement with the most-favored theory discussed 
above. As an example of the latter case, Chandra (1963) finds that the attachment- 
loss and diffusion terms are of equal importance near the F2 peak by a computation 
involving n e and q from Figure 8 and a loss coefficient behavior given by Ratcliffle, 
Schmerling, Setty and Thomas (1956). 



Other empirical models have been developed to describe the electron density 
distribution of the F region from a knowledge of the electron density and an arbi- 
trarily-defined scale height both measured at the F2 peak. These models vary in 
complexity from the use of a simple Chapman distribution with a constant scale height 
(Wright, 1960) to a variable scale height (Yonezawa and Takahashi, 1960) and 
finally a variable scale height gradient (Chandra, 1963). The latter model, applicable 
to an isothermal region where 0 + dominates, provides good empirical agreement with 
all available mid-latitude profiles attempted at height ranges between 100 km below 
the peak to about 700 km. Chandra’s empirical model is illustrated in Figure 9 where 
the abcissa (N/N 0 ) is the ratio of the electron density at a given altitude to its value at 
the F2 peak, and the ordinate z is given by 

z 9 - 0 h~h m )/H 0 , (7) 

where (h-h m ) is the height measured from the F2 peak and H 0 is the scale height 
corresponding to atomic oxygen. The parameter a is a measure of the departure from a 
simple Chapman distribution and is given by 

« - (H 0 -H m )/H 0 , (8) 

where H m is the “scale height” at the F2 peak. Chandra finds that the “scale height” 
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values computed from available electron density profiles correspond closely to neutral 
gas scale heights taken from the 1961 COSPAR International Reference Atmosphere. 

7. The Helium Ion Layer and the Protonosphere 

For a given value at a reference altitude, the upper ionosphere electron density which 
is given by Equation (6) will be governed by the mean ionic mass and the sum of the 
electron and ion temperatures. Charged-particle density profiles obtained at mid- 
latitudes under quiet solar conditions predominantly show that T e + If is practically 
constant at altitudes above 450 km. Consequently, changes in scale height of the 
electron-ion gas in the upper ionosphere predominantly reflect the transition from one 
ion species to a lighter one. It will be shown in this section that early concepts of an 
upper ionosphere characterized only by 0 + and H + ions and by a reasonably constant 
transition altitude now should be discarded in favor of a three-constituent medium 
having an intervening helium ion layer with a thickness which varies considerably 
with the atmospheric temperature. 

Drag observations on the Echo satellite led Nicolet (1961) to first suggest that 
helium should be considered an important constituent of the upper ionosphere. The 
first experimental evidences for the existence of helium ions in the upper ionosphere 
were obtained independently by an indirect method (Hanson, 1962a) and by direct 
measurement (Bourdeau, Donley, Whipple and Bauer, 1962). Hanson, working 
from an ion density profile obtained by Hale (1961) concluded that the thickness of 
the helium ion layer is of the order of 2000 km, extending from 1200 to about 3400 km, 
and that the measured scale height for helium ions corresponds to an atmospheric 
temperature of 1600° K. From the observed helium ion content (10 12 ions/cm 2 ) and an 
assumed rate of ionization, Hanson computed an equilibrium time for helium ions of 
about three days. He thus postulated that there should be no large diurnal variations 



Fig. 10. Direct detection of helium ions from the Explorer VIII retarding potential experiment. 
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in the helium ion layer. This conclusion has not been corroborated by other experi- 
mental observations and by the theoretical work of Bauer (1962, 1963). 

The Explorer VIII results (Bourdeau et at ., 1962) relative to direct detection of 
helium ions are summarized in Figure 10. This retarding potential experiment differs 
only in geometry from the device used by USSR investigators (Krassovsky, 1959) on 
Sputnik III. The sensor consists of three electrodes arranged in planar concentric 
geometry. The inner grid is biased negatively so that the collector is responsive only to 
the flow of positive ions from the ionosphere. Because the satellite velocity largely 
exceeds the thermal velocity of the ions, the latter have a kinetic energy relative to the 
spacecraft in proportion to their mass. In accordance with an expression given by 
Whipple (1959) the resulting volt-ampere curve, illustrated in Figure 10, will be 
characterized by an inflection point for a mixture of 0+ and He+ and by a plateau for a 
mixture of 0 + and H + . The abscissa of the volt-ampere curve is the total retarding 
potential given by the sum of the satellite and applied collector (</> c ) potentials while 
the ordinate is the ratio of the ion current at a given retarding potential to its value at 
zero potential. Planar ion traps have a good history for positive ion density determina- 
tion but generally yield too high an ion temperature (Hanson and McKibbin, 1961 ; 
Bourdeau et al , 9 1961), possibly due to an enhancement in the normal component of 
the particle energy caused by non-planar sheath geometry (Hinteregger, private 
communication). This factor is not considered in Whipple’s theoretical relationship. 
However, because of the relative insensitivity of the shape of the volt-ampere 
curve to ion temperature, the ability to determine ion compositon is not altered 
significantly. For daytime conditions, the Explorer VIII observations during the 
active life of the satellite (November-December, 1960) were that 0 + predominated 
below 1100 km and that He+ was the principal ion from 1100 km to the maximum 
altitude (1600 km) of the measurements. Unfavorable vehicle orientation prevented 
valid measurements for nighttime conditions and at altitudes above 1600 km. 

Subsequent to the reporting of the first experimental evidence for the importance of 
upper atmosphere helium ions, Bauer (1962) theoretically showed that the thickness 
of the layer would change with atmospheric temperature even for the simple case when 
the relative concentrations are invariant at a reference altitude. Assuming that the 
ionic species are distributed according to diffusive equilibrium (Dungey, 1955; 
Mange, 1960), the geopotential altitude ( fi'y ) at which two ionic species have equal 
concentrations has been given by Bauer (1962) as 


h'ij — h ' 0 = Htj In m 3 . 

titude at which the coi 
Hi] kT / (mt — Wj)g 0 , 


where bt Q is the geopotential altitude at which the concentrations of these ions are in 
the ratio mj, and where 


( 9 ) 

: in 

( 10 ) 


T being the atmospheric temperature and g 0 the acceleration of gravity at the earth’s 
surface. 

Bauer (1963) has proposed a more detailed theoretical model for the thickness of 
the helium ion layer as a function of atmospheric temperature by considering the case 
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where the relative concentrations of the three principal constituents vary at the 
reference altitude. He points out that the relative concentrations at a given altitude 
should vary with temperature in the same proportion as the ratio of the corresponding 
neutral constituents. Using the temperature dependent escape rate given for hydrogen 
(Bates and Patterson, 1961) and for both hydrogen and helium (Nicolet and 
Kockarts, 1962) and Equations (9) and (10), he computed two transition altitudes as 
a function of atmospheric temperature, a lower transition altitudes where 0 + and He + 
ions have equal concentrations and an upper transition altitude where He + and H + are 



Fig. 1 1 . Comparison of experimental measurement with theoretically-obtained model (Bauer, 1 963) 
of ion transition altitude as a function of atmospheric temperature. 

of the same density. These curves are shown to be in reasonably good agreement with 
experimental observations in Figure 11. The circles in Figure 11 refer to observed 
altitudes at which «(He + ) = w(0+) while the squares represent measured altitudes at 
which w(He + ) = n(H + ). The corresponding temperatures were obtained on the same 
spacecraft and represent either direct measurements of T e or 7* or values of ( T e + 7*)/ 2 
computed from the observed electron-ion scale height. The abscissa, then, corresponds 
to the neutral gas temperature only if it is assumed that the electron, ion and neutral 
gas temperatures all are in equilibrium at the altitudes indicated. 

Experimental points 1 and 2 are taken from Bourdeau et al. (1962) and Hanson 
(1962a), respectively. Experimental points 4 represent the first preliminary results 
from the Ariel satellite (Willmore, Boyd and Bowen, 1962). The Ariel investigators 
found that under daytime conditions helium ions predominated at altitudes between 
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950 km and to at least 1200 km, the satellite apogee. For a nighttime condition, it was 
found that the base of the helium ion layer was below 650 km and that hydrogen ions 
dominated at 1200 km. 

The Ariel experiment is based on the same principle as that of the Sputnik III and 
Explorer VIII retarding potential experiments and yet, in the opinion of this author, 
represents a major advance in the development of ion spectrometers for upper iono- 
spheric studies. Because of the spherical geometry used, it should not be as sensitive to 
the direction of the ion trajectory and thus provides valid measurements of ion 
temperature in contrast to the planar ion traps. Additionally, the Ariel investigators 
set up the experiment so that the second derivative of the volt-ampere curve is tele- 
metered directly. This reduces the effects of slowly-varying vehicle potentials and 
other ion sheath parameters. These two basic differences make for a higher ion com- 
position resolution than is possible with the Sputnik III and Explorer VIII ion traps. 

The Ariel satellite observations are extremely important in explaining the failure 
(Ulwick and Pfister, 1962) to detect He + ions from scale-height changes in nighttime 
electron density profiles. As Bauer correctly points out, the thickness of the helium 
ion layer is smaller than a helium ion height at low temperatures characteristic of a 
nighttime profile. As a result, the detection by indirect methods requires an extremely 
accurate charged-particle density profile. The detection of this layer is expected to 
become even more difficult as we approach the year of minimum solar activity. 

Experimental points 3 and 5 in Figure 1 1 were obtained by the indirect method 
from electron and ion density profiles obtained on Scout rockets. These data are 
presented in the succeeding section to illustrate the diurnal variation of electron 
density in the extreme upper ionosphere. 

8. The Diurnal Variation of Electron Density in the Upper Ionosphere 

The diurnal variation of electron density in the upper ionosphere is a complicated 
function of the electron density at and altitude of the F2 peak, the atmospheric 
temperature, and the extremely variable heights of the two ion transition altitudes. 
This variability can be illustrated by comparison of two charged particle density 
profiles taken near the diurnal maximum and minimum during this portion of the solar 
cycle. That representative of diurnal maximum conditions is illustrated in Figure 12. 
The circled points are electron density values obtained by Bauer and Jackson (1962) 
from Faraday rotation observations at 72 Me. The left-hand ordinate is geopotential 
altitude which takes into account the altitude dependence of the acceleration of 
gravity while the right-hand ordinate is the true or geometric altitude. Hanson 
(1962a) and Bauer (1962) have developed expressions based on an expansion of 
Equation (6) for the electron density distribution in a multi-constituent isothermal 
medium in diffusive equilibrium. These expressions suppose a given ratio of the con- 
centrations of the pertinent ions at a reference altitude. In Figure 12, two such 
theoretical curves (Bauer, 1962) are shown to illustrate the dependence of the 
electron density distribution on the ionic constituents. It can be seen that the circled 
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points are more consistent with a theoretical curve for a mixture of 0 + and He + (solid 
line) than a mixture of 0 + and He + (dashed curve). 

The charged particle density profile obtained near diurnal minimum conditions is 
illustrated in Figure 13. Both ion density (Donley, 1963) and electron density were 
measured. In the latter case, data were obtained only to a geopotential altitude of 



Fig. 1 2. Indirect detection of helium ions from electron density profile obtained on a Scout rocket. 

700 km. The good agreement of the experimental n+ values with the electron density 
obtained by rocket measurement and with the value of n e measured by an ionosonde 
at the F2 peak speaks well for the validity of the data. 

It can be shown from the theoretical curve provided by Bauer (curve A in Fig- 
ure 5) that the experimental ion density profile corresponds most closely to a three- 



Fig. 13. Scout rocket measurements of the diurnal variation of the upper ionosphere. 
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constituent mixture rather than 0 + and H + . The closest lit of the experimental points 
to the theoretical curve in the ion transition altitude leads to an estimate of less than 
300 km for the thickness of the helium ion layer. The discontinuity in the experimental 
n+ values at geopotential altitudes between 900 and 1100 km probably can be ex- 
plained from the observation that the spacecraft went from a dark to a sunlit condition 



Fig. 14. Qualitative model of the detailed diurnal variation of the upper ionosphere. 

in this region. The discontinuity probably represents a transient departure from 
temperature equilibrium. 

Curve B in Figure 13 is the electron density profile replotted from Figure 12 to 
illustrate an observed diurnal variation of charged particle density in the upper iono- 
sphere. The observation that the nighttime electron densities actually exceed those for 
daytime conditions at geometric altitudes greater than 1000 km should not be sur- 
prising since the lower ion transition altitude and the large scale height for protons 
overcompensate for the higher daytime atmospheric temperature and the higher day- 
time electron density at the F2 peak. 

Sayers, Rothwell and Wager (1962) claim to have discovered a new ionization 
ledge with peak densities located at geometric altitudes varying from 1200 km at mid- 
night down to 600 km at midday. The apparent negative scale heights gradients which 
were reported carry an implication of an ionization source in this altitude region. 
Discontinuities in charged particle density profiles perhaps can be expected at sunrise 
and at the equatorial anomaly, yet the discovery of this ledge is claimed throughout 
the day and for all latitudes traversed by the Ariel satellite. Sayers et al. constructed 
their profiles from rf probe measurements of electron density at the satellite position. 
The nature of the Ariel orbit during the time the data were obtained was such that the 
altitude profiles are not time independent but rather represent observations taken over 
a change of several hours in local mean time. An alternative explanation, then, is that 
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these ledges are only apparent in that they do not represent altitude discontinuities but 
rather changes in ionization brought about by the difference in local time at various 
points in the Ariel orbit. This alternative explanation is illustrated qualitatively in 
Figure 14. The experimental curves in Figure 13 already illustrate the diurnal variation 
of the electron density in the upper ionosphere. The four curves shown in Figure 14 are 
inferred profiles for the hours indicated. They were constructed from a typical varia- 
tion of the electron density at the F2 peak as observed for May (the month of the Ariel 
observations) by ionosondes, an atmospheric temperature variation taken from 
reference atmospheres, and the ion transition model of Bauer illustrated in Figure 11. 
It can be shown for a satellite which requires several hours of local time to traverse the 
400-1200 km altitude region, that an apparent ionization ledge would be observed 
when the satellite passes through the higher altitude region at night and conversely 
when it passes through the lower altitude region during midday. Parenthetically, an 
interesting feature to note from Figure 14 is the approximate constancy with diurnal 
time of the electron density at a geopotential altitude of about 800 km (approximately 
1000 km geometric). 

9. The Altitude Dependence of Charged Particle Temperature 

It is of considerable importance to compare electron (F e ), ion (7*) and neutral gas 
( T g ) temperatures as a function of altitude because their interdependence perhaps is 
the most sensitive index of complex reactions, the understanding of which is the 
ultimate goal of upper atmosphere physicists. For the time being, it will be assumed 
that Ti and T g are equal. 

Because direct and indirect measurements of charged particle temperatures have 
been made under radically different conditions and because of the limitations of the 
kinetic gas temperature models, various investigators have provided conflicting 
answers to the important question of temperature equilibrium between electrons and 
heavy constituents under daytime conditions. It generally is accepted that temperature 
equilibrium predominates at all altitudes during the night under quiet solar conditions. 
In this section it will be shown for a restricted set of latitude and atmospheric con- 
ditions that during midday temperature equilibrium exists below 150 and above 450 
km, in accordance with recent reviews (Bauer and Bourdeau, 1962; Bourdeau and 
Bauer, 1 962). It will be shown, on the other hand, that there is a substantial difference 
between spaceflight results and those from some ground-based radar incoherent back- 
scatter experiments. 

Before proceeding, it is important to define the term “temperature equilibrium”. 
Actually because in the ionization process the electrons are created with high initial 
energies, their temperature will always be somewhat higher and only can approach T g 
asymptotically in time depending on the efficiency of the energy transfer mechanism. 
Temperature equilibrium in this paper is defined as a difference between T e and T g 
which is smaller than the uncertainties in reference atmospheres and in spaceflight 
experimental methods of measuring charged particle temperatures. It is estimated, 
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perhaps optimistically, that these uncertainties together are about ten percent of the 
absolute value of T g . 

Our knowledge of the behavior of the neutral gas temperature is based principally 
on atmospheric densities computed from satellite drag observations and an assumed 
atmospheric composition. Values of T g so derived must be considered approximate 



Fig. 15. Comparison of the altitude dependence of neutral gas and experimental and theoretical 
electron temperatures at mid-latitudes under quiet solar conditions. 


principally because of the observational gap in atmospheric composition. Plotted as a 
solid line in Figure 15 is the altitude dependence of the neutral gas temperature 
published by the U.S. Committee on Extension to the Standard Atmosphere 
(Sissenwine, 1962). The COES A reference atmosphere was selected because it depicts 
“typical mid-latitude year-round conditions averaged for daylight hours and for the 
range of solar activity that occurs between sunspot minimum and sunspot maximum”, 
conditions which most closely approximate those under which the charged particle 
temperature measurements have been made. 

The earliest theoretical study of the temperature difference between T e and T g was 
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advanced by Hanson and Johnson (1961). Because of improved cross-sectional and 
solar radiation data and because of the omission of an important process (electron 
impact excitation of vibrational levels of molecular nitrogen), this work now must be 
considered superseded by more recent theoretical models (Hanson, 1962b; Dal- 
garno, McElroy and Moffett, 1962). These models, which consider ionization by 
solar radiation only, involve several expansions (depending upon the energy transfer 
mechanisms) of the equation: 

2 q'E 

(7WV )/ t = J— (11) 

jflefc 

where q is the rate at which photoelectrons of energy E are released and r is the time 
which it takes for these electrons to transfer their excess energy in various processes. 
Considerations are given to loss of this excess energy by (a) inelastic ionizing collisions 
with neutral particles, (b) inelastic exciting collisions with neutral particles where the 
energy is stored either in metastable or vibrational states, (c) energy transfer to 
ambient ions and (d) elastic collisions with ambient electrons. In all but the latter case, 
much of the excess energy of the newly-created photoelectrons is removed before they 
selectively heat the ambient electron gas. 

Hanson (1962b) computed four models for the case where the daytime temperature 
in the isothermal region is 1200° K. These four models are based on two different 
altitude-dependent values for the electron production rate and two different assump- 
tions for the excess photoelectron energy. The two extremes of these four curves are 
plotted in Figure 15. It is seen that the maximum predicted departure from tempera- 
ture equilibrium proposed by Hanson corresponds closely to the altitude of maximum 
electron production illustrated in Figure 8. Possible small errors have been introduced 
in Figure 15 by assuming that the value of T e -T g actually computed by Hanson are 
invariant when T g in the isothermal region is raised from 1200° K as used in his 
reference atmosphere and the 1500° K value used here. 

At present, an experimental determination of the complete altitude dependence of 
T e -T g can be accomplished only empirically by comparing Langmuir probe measure- 
ments of electron temperature obtained on vertical sounding rockets with a reference 
atmosphere. Langmuir probes have required considerable development in order to 
overcome problems associated with the disturbance introduced into the medium by a 
conducting body. Extremely high electron temperatures were reported from observa- 
tions of large potentials on the Sputnik III spacecraft (Krassovsky, 1959). However, 
it is dangerous to infer electron temperatures by this indirect method because of the 
probability of a vehicle potential which is artificially controlled by differing surface 
work functions, for example, than by ambient ionospheric parameters. 

It wasn’t until 1961 that electron temperatures close to accepted kinetic gas values 
were first reported for the E region by Japanese investigators (Aono et al 1961) and 
for the upper ionosphere by use of the Explorer VIII Satellite (Serbu, Bourdeau and 
Donley, 1961). However, these data by themselves do not provide a complete altitude 
profile of T e . The Japanese data, which represent an average from two separate 
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devices flown simultaneously, were obtained over the altitude interval between 100 
and 170 km and are plotted in Figure 1 5. A separate rocket (NASA 6.04) containing an 
ejectable symmetric bipolar probe was launched at approximately the same latitude 
and time. These results (Brace, 1962) are plotted as the dotted curve in Figure 15. 
Parenthetically, they represent a significant modification by Brace of the data which 
was first reported in a preliminary form (Spencer, Brace and Carignan, 1962). The 
error spread of the overall curve has been estimated by Brace to be ± 10 percent. 
Also plotted in Figure 15 is the atmospheric temperature computed on the assumption 
of temperature equilibrium by Jackson and Bauer (1961) from their experimental 
electron density profile. This value, applicable to the region above 450 km, was 
selected because it was taken under identical solar conditions (indexed by observed 
solar decimeter flux observations) as the two electron temperature profiles. Consider- 
ing the uncertainty in the reference atmosphere and the experimental error flag, it can 
be concluded for these particular solar conditions that temperature equilibrium exists 
below 150 km and above 450 km. Also there is reasonable agreement between the 
experimental results and the theoretical work of Hanson. 

Dalgarno et al. (1962) have proposed four theoretical electron temperature 
models for a given noontime value of T g . The end result in the form of the altitude 
dependence of T e -T g for T g = 2000° K is presented in Figure 16. The alternative 
curves are based on some reduction of excess photoelectron energy by (a) ionizing 
collisions alone, (b) including possible energy transfer through exciting collisions 
leading to the metastable state, (c) excluding (b) but including possible energy transfer 
through exciting collisions leading to the vibrational state and (d) a combination of all 
energy transfer mechanisms. Also plotted in Figure 16 are experimental values for 
T e -T g taken from Figure 15. It is tempting from the agreement of the experimental 
data with curves a and b in Figure 16 to suggest that excess photoelectron energy in 
the lower altitude regions is removed principally by a combination of ionizing colli- 
sions and exciting collisions leading to the metastable state. 

Hanson (1962b) has advanced a theoretical possibility based on low thermal 
contact between the ions and neutral particles that the ion temperature at very high 
altitudes could depart from the neutral gas temperature and approach the electron 
temperature. He estimates that such a transition could possibly occur in the 600- 
1000 km region. It is extremely difficult to corroborate this possibility from experi- 
mental observations because the differences between T e , Tt and T g are small at these 
altitudes. The daytime electron density profile illustrated in Figure 12 shows that the 
sum of the electron and ion temperatures taken in both the He + and 0 + regions are 
approximately the same. The constancy of this sum with altitude also can be inferred 
in published data (King, 1963) from the Alouette Topside Sounder Satellite. This is 
evidence either that the differences in the electron, ion and neutral gas temperatures 
are vanishingly small above 450 km or that the ion temperature transition proposed 
by Hanson occurs above 1000 km, if at all. 

The agreement of Langmuir probe results with the theoretical work of Hanson 
and Dalgarno et al. suggests that, at midday, temperature equilibrium exists 
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below 150 km and above 450 km whenever ionization is the result of solar radia- 
tion alone. However , because of the small number of observations , generalization of 
this conclusion may be premature . Spencer, Brace and Carignan (1962) have made 
rocket flights at auroral latitudes and into a disturbed mid-latitude ionosphere. For 
these conditions they report E and F2 region electron temperatures which are syste- 
matically higher than the quiet mid-latitude results. The difference is perhaps re- 



Fig. 16. Theoretical (Dalgarno, 1962) and experimental (Brace, 1962) models of the altitude 
dependence of the difference between electron and neutral gas temperatures at mid-latitudes under 

quiet solar conditions. 

presentative of the influence of corpuscular radiation. There is also a small possi- 
bility that the electron temperatures reported for disturbed conditions by Spencer 
et al. (1962) are not representative of the most abundant electrons. Their symmetric 
Langmuir probe samples only those electrons which can overcome a retarding po- 
tential of a few tenths of a volt. Since a Maxwellian electron energy distribution has 
been established for quiet solar conditions, the data reported by Brace and plotted 
in Figure 15 and 16 should be representative of the entire electron population. 
However, the existence of a Maxwellian energy distribution under disturbed iono- 
spheric conditions has not been verified. 

The spaceflight observations in a quiet daytime ionosphere are consistent with the 
conclusions (Bowles, Ochs and Green, 1962) from ground-based radar incoherent 
backscatter experiments that the electron and ion temperatures are in equilibrium in 
the upper ionosphere except at sunrise and during disturbed conditions. Evans (1962), 
on the other hand, also using backscatter experiments reports different results. The 
principal differences for which no explanation can be given here are: (a) Evans 
reports temperature equilibrium at 200 km, an altitude at which both theory and the 
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spaceflight observations show a significant departure from equilibrium ; (b) the ratio 
Te/Ti is 1.6 and constant above 300 km as reported by Evans whereas both theory and 
the spaceflight observations show this ratio decreasing with altitude above 250 km ; 
(c) Evans reports that the sum of T e and T% increases with altitude up to 700 km 
whereas it is concluded from many rocket measurements of charged particle density 
profiles that this sum is constant with altitude in the upper ionosphere. The results 
obtained by Evans also are inconsistent with the diurnal variation of charged particle 
temperatures at altitudes above 450 km, which is discussed in the succeeding section. 

10. Diurnal Variation of Upper Ionosphere Temperatures 

In this section, a comparison will be made of the diurnal variations of charged 
particle and neutral gas temperatures for the isothermal altitude region above 450 km. 
At these altitudes our knowledge of T g generally is derived from satellite drag ob- 
servations and an assumed atmospheric composition. The drag observations show 
that density variations are correlated with solar activity. Although not the source 
of upper atmosphere heating, solar decimeter radiation which is observable at the 
earth’s surface is an indicator of this interrelationship. For a given level of solar 
activity Jacchia (1962), Harris and Priester (1962) and Paetzold (1962) conclude: 

(1) The neutral gas temperature has its maximum value in mid-afternoon; 

(2) This maximum would occur at sunset if absorption of solar ultraviolet radia- 
tion was the only source of atmosphere heating; 

(3) Consequently, there is a second source of heating associated with the solar 
wind; 

(4) These conclusions apply to the atmosphere at low latitudes. 

In the preceding section it was shown for a specific set of observations taken at 
mid-latitudes that the electron and neutral gas temperatures have approximately the 
same value above 450 km, which is consistent with theoretical computations based on 
solar ultraviolet radiation as the sole ionizing source. If the low latitude ionization 
generally can be explained entirely by solar ultraviolet radiation and if, as claimed, 
the diurnal variation characteristic of the neutral atmosphere requires a second heat 
source, then the additional heating mechanism should be one which is least likely to 
selectively heat either the light or heavy atmospheric constituents. The heating by 
hydromagnetic waves proposed by Dessler (1959) probably meets these criteria. It is 
not possible from the charged particle temperature data presented in this section to 
infer a diurnal neutral gas temperature variation precisely enough to confirm or deny 
the hypothesis of a second heat source which is effective at mid-latitudes. 

There is some variance in reference atmospheres in the absolute magnitude of the 
temperature at diurnal maximum and minimum. The model attributed to Priester by 
Jastrow (1961) which is plotted in Figure 17 has an average daytime value consistent 
with the COESA standard atmosphere. It corresponds to a level of solar activity for 
which the observed solar decimeter (10.7 cm) radiation flux is 150 w/m 2 /cps. Other 
models (Harris and Priester, 1962; Jacchia, 1962) generally show somewhat lower 
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temperatures at diurnal maximum. Also plotted in Figure 17 are experimental values 
all taken at mid-latitudes of (a) the neutral gas temperature measured directly by 
Blamont (1961, 1962) using rockets containing sodium vapor release experiments, 
(b) direct measurements of electron temperature from the Explorer VIII Langmuir 
probe (Serbu et al., 1961) and (c) values of T g computed from rocket measurements of 
electron scale height on the assumption of temperature equilibrium, for which 
Tg = i (T e + Ti). These data were obtained at various times and consequently have 



Fig. 17. Diurnal variation of upper atmosphere temperature for altitudes above 450 km at mid- 
latitudes, solar decimeter flux — 150 x 10 -22 w/m 2 /cps. 


been normalized by this author to the same solar activity (10.7 cm flux — 150 w/m 2 /- 
cps) using the empirical relation given by Harris and Priester (1962). The temper- 
atures derived from electron or ion profiles are taken (1) from Donley (1963), (2), 
(8) and (9) from Berning (private communication), (3) from Bauer and Jackson 
( 1962), (4) from Jackson and Bauer (1961), (5) from Hanson (1962), (6) and (10) 
from Knecht and Russell (1962), and (7) from Hanson and McKibbin (1961). 

It is tempting from the general agreement of the neutral gas temperatures with the 
Explorer VIII electron temperatures and with the values of \{T e + Ti) derived from 
charged particle density profiles, all illustrated in Figure 17, to infer support for the 
conclusion reached in the previous section that temperature equilibrium exists above 
450 km for most of the day under quiet ionospheric conditions at mid-latitudes. How- 
ever, the empirical method of comparison used in Figure 1 7 can only be an approximate 
test of equilibrium, because except for Blamont’s data, the neutral gas temperature is 
an inferred parameter and because the preciseness of normalizing data taken at 
different times to the same level of solar activity has not been established. 
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The Explorer VIII electron temperature device is a gridded trap which experi- 
mentally eliminates unwanted photocurrent effects. The error flags shownin Figure 17 
result from a limited volt-ampere curve resolution imposed by the telemetry system. 
There possibly is an additional small error associated with the electrical transparency 
characteristics of the grids used to remove the photocurrent effect (Serbu et aL, 1961). 
The absence of an on-board tape recorder prevented detailed analyses of the diurnal 
electron temperature variation and the detection of latitude effects. 



Fig. 18. Possible latitude dependence of the diurnal variation of electron temperature for altitudes 
above 450 km, solar decimeter flux m 100 X 10~ 22 w/m 2 /cps. 

A different type of Langmuir probe now has been flown on the Ariel satellite. 
Basically an asymmetric probe with the sensor located on a boom away from the main 
spacecraft, the experiment tends to suppress such unwanted effects as varying vehicle 
potentials and photocurrent and magnetic field modulation of the electron current by 
measuring the electron temperature in an extremely short time interval. The avail- 
ability of an on-board tape recorder has permitted the investigators to report the most 
detailed diurnal electron temperature variation yet obtained and additionally to 
separate the data according to latitude. Their preliminary results as reported (Will- 
more, Boyd and Bowen, 1962) are shown in Figure 18. The three curves, each rep- 
resenting a different latitude region were obtained at a time when the 10.7 cm solar 
radiation flux averaged 100 w/m 2 /cps. The Priester model of the diurnal variation of 
the neutral gas temperature in Figure 17 has been normalized to this lower level of 
solar activity and replotted in Figure 18 for comparison. 

The Ariel data were analyzed on the assumptions that the temperature above 
450 km is isothermal and that there is a symmetrical distribution about the equator. 
Willmore et aL point out that the salient features are a general rise in temperature 
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and a less pronounced diurnal variation at high latitudes and that both of these 
observations are consistent with significant atmospheric heating by particles being 
dumped at high latitudes. 

Another salient feature of the Ariel data is the high electron temperature observed 
in the sunrise period. This is not inconsistent with the experimental results in Figure 17 
because of the observational gap in the latter data during this period. The sunrise effect 
in the Ariel data is evidence for a departure from temperature equilibrium. The diurnal 
shapes of the curves suggest that temperature equilibrium is re-established toward 
midday. The observation of a departure from temperature equilibrium at sunrise is 
consistent with the conclusions from one set of ground-based backscatter results 
(Bowles, 1962) but not with Evans (1962) who reports maximum departure from 
equilibrium at noon. A higher ratio TejTi at sunrise than at noon is more consistent 
with what would be expected from a cursory extension to the theoretical observation 
(Hanson, 1962) that T e /Ti is directly proportional to qjn e 2 above the F2 peak. The 
electron production rate, g, is in turn proportional to atmospheric density. It can 
be shown from a comparison of the atmospheric density variation derived from 
satellite drag with the electron density profiles illustrated in Figure 14 that the ratio 
q/n e 2 and thus that T e jTi should be higher at dawn than at noon. 

It would be premature to draw firm conclusions regarding additional ionizing and 
heating sources from the diurnal comparisons of upper ionosphere electron and 
neutral gas temperatures shown in Figures 17 and 18. The experimental charged 
particle observations in Figure 17, for example, were obtained over a period of more 
than one year. The accuracy of the process of normalizing these data to the same level 
of solar activity quite likely introduces enough error that the diurnal shape of the 
indirectly-derived neutral gas temperature cannot be confirmed. Consequently, the 
hypothesis of a second heat source at low latitudes also cannot be confirmed. 

The midday values of T g proposed by Priester and which is applicable to low 
latitudes fall in between the two lower latitude electron temperature curves measured 
on the Ariel satellite. No conclusion can be drawn from the different diurnal shapes of 
T e and T g in Figure 18, principally because two refinements in the analysis of the Ariel 
data are required. Firstly, there is no evidence that the Ariel results have been nor- 
malized for the day-to-day fluctuations in solar activity which can amount to a ten 
percent correction. This normalization could alter both the diurnal shape of the 
individual electron temperature curves and the amplitude of the latitude effect. 
Secondly, it would be expected again from a cursory theoretical analysis that the 
electron temperature in the upper ionosphere during the sunrise period would be 
altitude dependent. Since the present analysis of the Ariel data assumes no altitude 
effect, this also could alter the diurnal shape of the individual electron temperature 
curves. With reference to the diurnal variation of T g , the Ariel data suggest a latitude 
effect which up until now has not been seriously considered in neutral gas models. 
This latitude effect in turn suggests that other heat and ionization sources in addition 
to direct solar radiation need to be considered in a description of even the normal 
auroral ionosphere. 
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11. Summary 

Spaceflight observations made in recent years are most consistent with the following 
general conclusions: 

(1) The principal ionizing agents for the D region under quiet solar conditions at 
an epoch corresponding to the middle of the solar cycle are cosmic radiation for the 
region below 70 km and Lyman a radiation at altitudes between 70 and 85 km. The 
negative ion abundance at all D region altitudes has been underestimated theoretically. 

(2) X-radiation (2-8 A) and energetic protons, which are respectively responsible 
for Sudden Ionospheric Disturbances and Polar Cap Absorption Events, lead to a 
D region electron density enhancement of about two orders of magnitude. 

(3) The principal ions produced in the lower ionosphere are N 2 + , C> 2 + and 0 + . The 
ions which exist in greatest abundance are NO+, C> 2 + and 0+. These ionization char- 
acteristics result from the combined influence of solar ultraviolet and X-radiation 
in the E region with ultraviolet radiation becoming the dominant source in the F 
region. 

(4) The average E region electron density decreases by about two orders of 
magnitude at night. Some types of sporadic-E ionization have a depth of less than one 
kilometer and a horizontal dimension greater than 70 km. Other types have a smaller 
horizontal dimensions. There is a first order increase of the electron density in a 
sporadic E-layer. 

(5) The formation of the F2 peak is best explained by the combined altitude 
dependence of the electron production rate, electron loss by an attachment like process 
and diffusion. 

,(6) The constant electron-ion scale height which is observed for altitude intervals 
of a few hundred kilometers above the F2 peak is evidence that the upper ionosphere 
is isothermal and in diffusive equilibrium. 

(7) The upper ionosphere is characterized by the presence of three ionic con- 
stituents (0 + , He + , H+) each predominating in a different altitude region. The 
thickness of the helium ion layer has a large diurnal variability. The control which 
these ionic constituents exercise over the electron density distribution is such that the 
electron density is generally higher at midday below 1000 km, shows a small diurnal 
variation near 1000 km, and is generally higher at altitudes between 1000 and 2000 km 
at night. 

(8) In an undisturbed mid-latitude ionosphere, it is generally accepted that 
temperature equilibrium can be expected at night. Significant departures from equi- 
librium can be expected at sunrise and during disturbed conditions. Toward midday 
under quiet mid-latitude conditions significant departures from equilibrium would 
be expected in the lower F region, but small differences between T e , Ti and T g are 
observed below 150 and above 450 km. 

(9) Except for the sunrise anomaly the diurnal maximum and minimum values of 
electron temperature at mid-latitudes approximate the diurnal extremes of the neutral 
gas temperature in the upper ionosphere. An observed latitude dependence of electron 
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temperature suggest that corpuscular radiation is an important heat source at high 
latitudes. 
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APPENDIX 


ON THE INTERACTION BETWEEN A SPACECRAFT 
AND AN IONIZED MEDIUM 


1. Introduction 

Much of the data discussed in the preceding sections were obtained by direct meas- 
urements techniques (plasma probes). The success of experiments of this type is 
dependent upon an evaluation of the effects of the interaction between the spacecraft 
and the ionized atmosphere immediately surrounding it. An estimate of the errors 
which this interaction could introduce into the results can best be made by the 
respective investigators. However, it should be useful to readers unfamiliar with this 
class of experiments to present general considerations of this interaction and estimates 
of the relative confidence that can be placed on the electron and ion parameters 
obtained from plasma probes. 

2. The Effects of RF Fields and of Contaminant Outgassing 

Two factors artificially introduced by the spacecraft are the presence of radio- 
frequency fields due to telemetry transmissions and outgassing phenomena. The 
Japanese work (Aono et ah , 1962) on resonance probes demonstrates that if the 
frequency of rf transmissions exceeds the plasma frequency, there is a negligible 
perturbing effect. Whale (1963) points out that this is so only for low values of 
radiated power. Low enough radiated rf power at a high enough frequency was 
maintained on the spacecraft from which data have been reported in the preceding 
sections of this paper to conclude that the rf field can be neglected. 

The agreement of a direct satellite measurement of atmospheric density (Sharp, 
Hanson and McKibbin, 1962) with that deduced indirectly from drag observations 
provides confidence that contaminants outgas soon enough to conclude a negligible 
effect from this factor on plasma probe results from orbiting vehicles. However, this 
effect must be taken into consideration for results obtained from vertical sounding 
rockets. Ulwick and Pfister (1961), for example, find that outgassing contaminants 
actually mask the desired results for the case where care is not taken to seal the 
rocket’s propellant tanks after engine burnout. Of the vertical sounding rocket data 
considered in this paper, that of Brace (1962) presented in Figures 15 and 16 are 
expected to be most free of errors due to outgassing phenomena because of the care 
taken in separating the Langmuir probe from the rocket booster. The Japanese 
appear to have overcome the problem by placing their probes on booms extended 
a considerable distance from the rocket body. Of the rocket data discussed in this 
paper, those of Smith (1961a, 1961b) shown in Figures 3 and 7 are most susceptible 
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to outgassing effects since he uses an active electrode flush-mounted on the rocket 
body. Even here, however, the good agreement which he obtains between ascent and 
descent data suggests a minimal error from the outgassing effect. 

The remaining disturbance factors can be summarized as the sensitivity to space 
craft orientation of : (a) the potential of the spacecraft relative to the ambient medium- 

(b) the current (/+) due to diffusion of positive ions from the medium to the spacecraft, 

(c) the current (i e ) due to diffusion of electrons from the medium to the spacecraft 
and (d) photocurrent (i p ) due to solar-radiation induced emission of electrons from 
the spacecraft surfaces. An experimental model of these interaction factors has been 
obtained from a set of experiments flown on the Explorer VIII satellite (Bourdeau 
et al ., 1961). This model is qualitatively consistent with that proposed by Gringauz 
and Zelikman (1957). The Explorer VIII results are summarized herewith to permit 
subsequent general consideration of possible errors in the determination of iono- 
spheric parameters by direct measurements techniques. 


3. The Potential of a Spacecraft in the Ionosphere 


The average equilibrium potential of a conducting body at rest where magnetic fields 
and solar radiation may be neglected is given by 



Je 

J+ ’ 


( 1 ) 


where J e and J+ are the respective random electron and positive ion current densities 
which flow in the plasma in the absence of the conducting body. The respective 
current densities (J e and »/+) are proportional to the charged particle density and the 
square root of the charged particle temperature and inversely proportional to the 
square root of the mass of the charged particles. Thus, for a condition near temperature 
equilibrium, the ratio J e /J+ is proportional to the square root of the ratio of the 
electron and ion masses. Because of the much smaller electron mass, J e > J+ and (f> 0 
will be negative for the stated conditions of a motionless body at night with no 
magnetic field. 

For a moving body at night with no magnetic field, the average equilibrium 
potential is given by 


</> o = — • 


kTe . Je Se 


In 


J+ S+ 


( 2 ) 


where S e and S+ are the areas over which the respective electron and ion currents 
are effective. Because of the higher electron mobility, S e > S+. Since J e already is 
much greater than /+, (f> 0 will be even more governed by the characteristics of the 
ambient electrons, specifically T e . Indeed, the observation of a <f> 0 of a few tenths of 
a volt negative for a nighttime condition on the Explorer VIII satellite is consistent 
with Equation 2 (using reasonable values of T e ) and thus lends confidence to a 
conclusion that the basic laws of kinetic theory are being obeyed for the case of an 
uncomplicated spacecraft. 
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For a daytime condition, the photocurrent effect which is of the same polarity 
as the current due to diffusion of positive ions from the medium to the spacecraft 
tends to drive the vehicle potential positive such that, for no magnetic field 


For a specific value of the ambient charged particle density, the diffusion and photo- 
current effects can be expected to cancel with the convenient end result that </> 0 is zero. 
This was observed on the Explorer VIII satellite for the case where N e = 10 4 cm -3 , 
a value consistent with measured values of J p and values of J e computed from 
separate measurements of N+ and T e . 

The preceding discussions of spacecraft potential were based on the simplifying 
assumption of no magnetic field. The motion of a spacecraft with velocity V through 
a magnetic field B will cause the vehicle surface to be non-equipotential. As a result, 
the vehicle potentials quoted above apply only to those measured on a plane passing 
through the center of the spacecraft and perpendicular to the vector representing the 
cross-product of the velocity and magnetic field vectors. The more general term for 
the potential of any surface element on a spacecraft is 

<f> — (f>o + (V xB).d, (4) 

where d is the vector distance from the surface element in question to the center of 
the spacecraft. A magnetic-field induced emf of 0.14 volts was measured across 30 cm 
on the Explorer VIII satellite (Bourdeau et al., 1961). This is consistent with a 
computation involving the insertion of known values of V, B and d into Equation (4). 
More recently, Willmore et al. (1962) have reported a magnetic-field induced emf 
of 0.5 volts measured across a larger distance on the Ariel satellite. 

The preceding discussion demonstrates that the behavior of the potential of a 
surface element on a spacecraft with a smooth configuration and without extensive 
protuberances is somewhat predictable. It can be expected that the potential of a 
surface element will undergo time changes which reflect different amounts of the 
overall vehicle area being exposed to various effects, particularly that due to the 
magnetic field. Thus as far as the effect of vehicle orientation modulation of spacecraft 
potential is concerned, it is necessary either to make a measurement during a time 
interval short compared to the orientation period or to resort to statistical averaging 
in extracting those geophysical parameters which are sensitive to the spacecraft 
potential. 

4. The Reliability of Ion Parameter Measurements 

In addition to the potential of the vehicle surface element at which the measurement 
is made, the other controlling factors in the determination of positive ion parameters 
is associated with velocity modulation of the ion current. Excepting the case of a 
measurement at the apogee of a vertical sounding rocket, the vehicle velocity greatly 
exceeds the thermal ion velocity. As a result, the current (i+) flowing to a surface 
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element of area A pointed within 45° of the velocity vector is approximately given by 

/+ = e A N+ V cos 9 , | 9 | < 45° , (5) 

where 9 is the angle between the surface element and the velocity vector. The larger 
vehicle velocity produces a “ram effect” such that i+ would be expected to be a 
maximum when 9 is zero and would vanish when the surface element is pointed into 
the wake of the vehicle. This was demonstrated on the Explorer VII l satellite by 
parametric scanning of /+ as the vehicle spun. The effect illustrated in the top graph 
of Figure 19 was measured by a gridded device where, for this case, unwanted electron 
and photocurrents were experimentally removed. Equation (5) is valid whenever the 
potential of the surface element at which the measurement is being made is negative 
with respect to the plasma. It is seen from Equation (5) that the positive ion density 
(7V+) is computed from known parameters rather than resorting to assumptions about 
other ionic properties. Accordingly, reasonable assurance can be had that whenever 
the conditions of a negative surface element and a low value of 9 are met, accurate 
values of N+ will result. The validity of ion density determination whenever these two 
conditions are met has been demonstrated in several ways. Firstly, there is good 
agreement (Donley, 1963) between a direct measurement of n+ and a separate but 
simultaneous radio-propagation method of measuring n e . This agreement is illustrated 
in Figure 13. Secondly, consistent conclusions relative to the ionic composition 
(cf. Figure 11) and charged particle temperatures (cf. Figure 17) are reached whether 
one uses directly-measured ion density profiles or electron density profiles obtained 
by radio-propagation methods. 

The history of the measurement of positive ion density by planar devices for the 
case of a negative surface element but where 90° < 9 < 45° or alternatively where 
the thermal ion velocity approaches or exceeds the spacecraft velocity is not as good. 
Whipple (1959) has set forth the working equations for a planar ion trap which takes 
into account all values of the spacecraft-to-ion velocity ratio: 

_ x 2 

/+ = eN+A [V cos 6 + \ erf x) + — ] , (6) 

2\/n 

where a+ is the ion thermal velocity and where 

V cos 9 cj> e 

x = 

a+ kT+ 

The above equations apply for the case where the ion sheath and the sensor used in 
making the measurement have the same geometry. There is experimental evidence 
that this assumption only holds approximately for the case of planar geometry. This 
was pointed out by Bourdeau et al (1961) from the observation that /+ for 9 = 90° 
is too large as illustrated in the top graph of Figure 19. Thus, ion density values 
computed for such angles from the planar ion trap on the Explorer VIII satellite 
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would be too high. Similarly, Brace (1962) finds that the values are ten percent 
too high at apogee of a sounding rocket although he gets good agreement with 
simultaneous measurements of electron density obtained by a radio-propagation 
experiment, at the lower altitudes where V > a+. 

Continuing the same reasoning, results show that Equation (6) cannot be applied 
exactly in determining other ionic parameters from planar sensors even when 0 = 0. 
For the condition of a ram ion current and a known vehicle velocity, it is possible 





Fig. 19. Experimental observations of current exchange between the ionosphere and the 

Explorer VIII satellite. 

to compute an idealized behavior for /+ as a function of the potential of a surface 
element and as the ionic parameters of mass and temperature are varied. Such an 
idealized behavior for a constant T+ with different values of m+ is shown in Figure 10 
where it is illustrated that a mixture of 0 + and He + is readily identifiable from one 
of 0+ and H + . It should be possible also to measure T+ from the slope of the curves 
in the regimes where the current behavior is due to a single ionic constituent. However, 
Equation (6) was derived by Whipple for the case of a planar sensor mounted on an 
infinite plane, a condition which only has been approximated on spacecraft flown to 
date. The exact conditions encountered are difficult to treat theoretically. It is most 
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likely because the geometry of the ion sheath is not planar that ion temperatures 
computed from planar ion traps by use of Equation (6) are generally too high. 
However, these considerations should not affect significantly the use of planar ion 
traps as ion composition measuring devices for the following reasons. The error in 
the measured ion temperature can be estimated from the difference between those 
experimentally-obtained by use of Equation (6) and estimates of the neutral gas 
temperature. It can be shown that within the limits of the error so obtained, the 
shapes of the curves illustrated in Figure 10 are not affected enough to alter the 
characteristic inflection point which identifies a mixture of 0 + and He + and the 
characteristic plateau which identifies a mixture of 0 + and HV The method for 
demonstrating this involves maintaining a constant ionic mixtiure and varying the 
ionic temperature within the limits of the estimated error in TV in Equation (6). 
It should be emphasized that the preceding criticism of planar ion traps for TV 
determination does not apply to spherical ion traps such as those used on the Sputnik 
III and Ariel satellites. To the contrary, it would be expected that an equation for 
spherical geometry analogous to Equation (6) could be used with confidence for all 
ionic properties because of the high probability that the sensor and its surrounding 
ion sheath do have the same geometry. This conclusion is supported by the fact that 
experimental values of 7V obtained by spherical devices are consistent with accepted 
values of the neutral gas temperature. These considerations in turn should permit the 
measurement of ion composition to a higher resolution than is possible with the 
planar ion traps. 

To summarize the validity of ionic properties, it is expected that the ion density 
profiles obtained by Hale (1961), Hanson and McKibbin (1961) and that reported 
by Donley (1963) in Figure 13 closely correspond to the true geophysical situations. 
This additionally is true of ion densities reported from Sputnik III. The ion composi- 
tion results represented in Figures 10 and 11 also are expected to be approximately 
geophysically correct. Until additional theoretical work has been done on the ion 
sheath problem, ion densities obtained either when the ion thermal velocity approaches 
that of a spacecraft or from planar traps pointed at too large an angle from the 
velocity vector are expected to be incorrect. This is similarly true of the ion temperature 
parameter when obtained from a planar trap irrespective of the pointing angle. None 
of the data obtained for the latter two conditions has been used to derive geophysical 
conclusions in this paper. Parenthetically, the complicated behavior of the potential 
at a particular vehicle surface element has been taken into account in summarizing 
those data considered valid. It can be seen from Figure 10 that for negative potentials, 
the ion current is independent of the potential. Since for all cases used in the preceding 
sections the potential was negative when ion density was computed, the values of the 
latter parameter should not be subject to errors due to a changing vehicle potential, 
including the magnetic field induced emf effect. In ascertaining ion composition from 
the Explorer VIII planar ion trap, one is concerned with the behavior of /+ as the 
potential is varied in a positive direction. Errors due to changes in overall vehicle 
potential were minimized, in this case, by obtaining a volt-ampere curve for changes 
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in orientation as little as 15 degrees (Bourdeau et al 1961). The same philosophy of 
deriving ionic parameters in a short time interval was used by Willmore et al. (1962) 
for the obtainment of Ariel satellite ion composition and temperature results. 

5. The Validity of Electron Parameter Measurement 

The derivation of electron parameters from plasma probes is a degree more difficult 
than the case of positive ions principally because of the higher sensitivity to spacecraft 
potential effects. Digressing from this point a moment, the electron current to a 
particular surface element of constant potential will be dependent on the vehicle 
velocity, although to a lesser extent than the ions. This is because the thermal electron 
velocity exceeds even the satellite velocity. Electron current as a function of orientation 
relative to the velocity vector was uniquely measured on the Explorer VIII satellite 
(Bourdeau et al 1961) by a device in which the positive ion and photocurrents were 
removed experimentally. These data are reproduced in the center graph of Figure 19. 
The amount of this modulation is consistent within a factor of two with what would 
be expected from an independent measurement of the electron temperature, the 
difference representing a gap in the theory of ion sheath behavior. This difference can 
be reflected in a measurement of electron density but should have a negligible effect 
on electron temperature determination. 

Before proceeding to a fuller discussion of the validity of electron parameters, it 
is necessary to illustrate the behavior of the total net current flowing to and from a 
spacecraft surface element as the vehicle rotates. The net current includes the effects 
of ambient electron and positive ion diffusion and photoemission. Such data obtained 
from the use of an insulated portion of the Explorer VIII satellite by Bourdeau et al. 
(1961) are reproduced in the bottom graph of Figure 19. Three major effects all of 
which mask velocity modulation of the electron current are to be noted. The first is 
the peak in the ram ion current as the sensor points in the direction of motion. The 
second and by far the major peak is due to photoemission as the sensor points toward 
the sun. The relative amplitude of this peak is dependent upon the ambient electron 
density which for this illustration was approximately 10 4 cm -3 . The third gross feature 
is the peak in the electron current which is due to the magnetic field effect. It occurs 
on the graph at a velocity azimuth angle between 270 and 360 degrees. At this position, 
the sensor was pointing in the V x B direction, the position of the most positive 
potential. For this reason, the peak in i e shown in Figure 19 is to be expected. 

All of the factors brought out by the graphs in Figure 20 must be borne in mind 
in the ensuing discussion of techniques for the measurement of electron parameters. 
The simplest Langmuir probe, yet the most difficult to analyze, is of the exposed 
asymmetric type. In the version used by Smith (1961a, 1961b), the collector is the 
rocket’s nose tip whose surface area is small enough that the rocket potential can be 
considered the reference potential. The observed collector current as a function of the 
collector-to-rocket potential ( <f > cr ) taken at an altitude of 115 km for a daytime 
condition is illustrated in Figure 20a. One needs to know the behavior of only the 
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incoming electron current to compute electron density and temperature. This requires 
an assumption as to the behavior of the unwanted ion and photocurrents to which the 
collector also is sensitive. The usual method is to assume that the current observed at 
high negative potentials is entirely due to i+ and i v and by extrapolation subtract 
their total behavior. The corrected curve illustrated in Figure 20b then is plotted 
logarithmically and electron temperature is computed from the slope in the exponen- 
tial regime. For this particular set of data, the reported electron temperature was 
1500° K, a value considerably higher than was obtained under quiet solar conditions 



(a) UNCORRECTED PROBE CURRENT (b) PROBE CURRENT CORRECTED FOR 

ELECTRON TEMPERATURE DETERMINATION 

Fig. 20. Typical volt-ampere curves (Smith, 1961a) from a flush-mounted asymmetric 

Langmuir probe. 


for the same altitude but by different devices (Aono et al . , 1961; Brace, 1962). The 
difference could be geophysical in nature since the data reported by Smith (1961a) 
were obtained within twenty-four hours of a geomagnetic disturbance. However, it 
also is possible that the variance is attributable to the different techniques used. 

It can be seen that the slope of the curve shown in Figure 20b is very sensitive 
to the assumed i+ and i v corrections. Arguments can be made that the effect on electron 
temperature determination still should be small because /+ and i v vary quite slowly 
with <f> cr whenever the collector is negative with respect to the plasma potential. The 
latter point is identified in Figure 20a as the value of <j) C r where the collector current 
has the value (/g) 0 . The assumption of a slowly varying /+ and i v is likely to be correct 
for probes extended away from a large rocket body. However, for the flush-mounted 
asymmetric probe used to obtain the data shown in Figure 20, there is a small pos- 
sibility of current exchange between the collector and the adjacent rocket surfaces. 
Evidence for such a possibility exists in the form of the shift of the overall curve from 
the origin, a bias which Smith attributes to dissimilar collector and rocket surfaces. 
It is possible, by reason of this bias, that some of the rapid change in collector current 
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near ( i e )o is due to unexpected large variations in the unwanted ion or photocurrents 
and that this could be reflected as an error in the electron temperature determination. 
The possibility of such errors can only be suggested, not verified at this time. 

It should be emphasized that these considerations do not affect significantly the 
electron density values reported by Smith and used in Figure 7 to illustrate the diurnal 
variation of the E region. Again referring to Figure 20, electron density is computed 
from the measured value of ( i e ) 0 according to 

(ie)o = (eN e c e A)/ 4, (7) 

where c e is the mean electron velocity which in turn is proportional to T e L It is seen 
from Figure 20 that the value of ( i e ) 0 is some ten times the observed total effect of the 
unwanted ion and photocurrents. It follows that the electron density determination 
should be insensitive to the methods used in correcting for /+ and i p even though the 
electron temperature determination is sensitive in this respect. Smith found that the 
electron density computed from Equation (7) by insertion of the electron temperature 
obtained from the exponential slope of the corrected curve is a factor of two higher 
than that obtained from an ionosonde. This factor could be explained if the slope of 
the corrected volt-ampere curve and thus the derived electron temperature were in 
error by a factor of 4. Alternatively, it could be explained by a stronger vehicle velocity 
effect than expected as was pointed out in the discussions accompanying Figure 19. 
In either event, since normalization to the ionosonde data did take place the only real 
error in question is the shape of the electron density profile. If, instead of the assump- 
tion of an isothermal electron temperature used by Smith in deriving the daytime 
profile shown in Figure 7, one were to assume an electron temperature behavior 
similar to that of the neutral gas then the errors in the profile would only be about 
20 percent. Similarly, the nighttime electron density results should be approximately 
correct. Consequently, although the electron temperature results obtained from the 
flush asymmetric probe are subject to some question, the important findings of Smith 
relative to the diurnal variation of E region electron densities are not significantly 
altered by these considerations. 

It should be emphasized that the suggestion that asymmetric probes flush-mounted 
on a rocket could be subject to errors in electron temperature do not necessarily apply 
to other Langmuir probes that have been flown on spacecraft. As far as the electron 
temperature parameter is concerned the resonance probe of Aono et al . (1961) is in 
principle an asymmetric Langmuir probe. However, its collector is extended away 
from the main rocket body which should suppress ion and/or photocurrent exchange 
between the collector and the main rocket body and additionally should make the 
results insensitive to perturbations of the ion sheath surrounding the main body. 
Similarly, confidence can be placed that these types of errors are suppressed in the 
ejectable asymmetric Langmuir probe used by Spencer et al. (1962) especially since 
care has been taken in providing appropriate guard rings. It must be noted, however, 
that the latter two devices measure electron temperature when the collector is negative 
with respect to plasma potential. For this reason, these probes could contain errors 
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if there existed the unlikely possibility that the electron energy distribution is not 
Maxwellian. The only other criticism that can be made of these devices is associated 
with the time taken to obtain the respective volt-ampere curves. When this time is * 

long, there is a possibility that the electron current is varying with vehicle velocity or 
magnetic field modulation (as shown in bottom graph of Figure 19) rather than 
entirely because of the application of the retarding potential. When the electron 
current varies with probe orientation rather than probe potential, errors in electron 
temperature can result. Because of the possibility of this particular type of effect, it 
is considered by this author that the significant changes in electron temperature over 
small altitude intervals reported by Spencer et al. (1961) are likely not to be real. 
However, the errors in the overall smoothed altitude curves which Spencer et al. 
produced by statistical averaging should be small, a conclusion applicable to the 
results of Brace (1962) used in Figures 15 and 16. 

The electron temperature data used in Figure 17 were reported by Serbu et al. 

(1961) who used a gridded device on the Explorer VIII satellite. By the use of grids, 
unwanted ion and photocurrents were separated out experimentally permitting un- 
ambiguous measurement of the wanted electron current. Additionally, volt-ampere 
curves were obtained for small changes in vehicle orientation. On the other hand, the 
use of grids introduces a source of second order errors - that due to changes in the 
electrical transparency of the grids with retarding potential. The device used by 
Willmore et al. (1962) to obtain the data shown in Figure 18 is an exposed probe t 

with appropriate guard ring extended on a boom. The transmission of the first and 
second derivative of the volt-ampere curves should make for electron temperature 
determination over extremely short intervals of time and for all pertinent values of 
probe-to-plasma potential. These features should permit the investigators to report 
data insensitive to errors due to orientation-induced changes in unwanted ion and 
photocurrents. 
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